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Abstract 



' Using the transport equations for an ideal anisotropic collisionless plasma derived from the Vlasov 

equation by the 16-moment method, we analyse the influence of pressure anisotropy exhibited by 



collisionless magnetized plasmas on the magnetothermal (MTI) and heat-flux-driven buoyancy 
(HBI) instabilities. We calculate the dispersion relation and the growth rates for these instabilities 
in the presence of a background heat flux and for configurations with static pressure anisotropy, 
finding that when the frequency at which heat conduction acts is much larger than any other 
, frequency in the system (i.e. weak magnetic field) the pressure anisotropy has no effect on the 



' MTI/HBI, provided the degree of anisotropy is small. In contrast, when this ordering of timescales 

does not apply the instability criteria depend on pressure anisotropy. Specifically, the growth time 
of the instabilities in the anisotropic case can be almost one order of magnitude smaller than 
its isotropic counterpart. We conclude that in plasmas where pressure anisotropy is present the 
MTI/HBI are modified. However, in environments with low magnetic fields and small anisotropy 
such as the ICM the results obtained from the 16-moment equations under the approximations 
considered are similar to those obtained from ideal MHD. 
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I. INTRODUCTION 



In many magnetized dilute astrophysical plasmas thermal electron conduction occurs 
almost exclusively parallel to magnetic field lines. In this regime, the equations of ideal 
magnetohydrodynamics (MHD) that describe the plasma dynamics must be supplemented 
with anisotropic transport terms for energy and momentum due to the near free-streaming 



motions of particles along magnetic field lines 



10| . Balbus {2,4] have shown that anisotropic 



thermal conduction can fundamentally alter the Schwarzschild convective instability crite- 
rion in stratified atmospheres, in that convection sets in when a temperature gradient (as op- 
posed to an entropy gradient) ispresent. This new convective instability is termed the mag- 



netothermal instabi 
Later on, Quataert 



ity (MTI) {2] when the temperature decreases with height {g ■ VT > 0). 



331 ] showed that in the presence of a background heat flux there appears 



a heat-flux-driven buoyancy instability (HBI) when the temperature increases with height 
{g ■ VT < 0). Essentially, the MTI and the HBI occur because the heat flux must follow the 
perturbed magnetic field lines, and thus there are regions in the plasma which are locally 
heated or cooled [see 



24 



26 



33|. 



The MTI and HBI couple the magnetic structure of the plasma to its thermal properties 



n 
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4l|. In a 



and can have important implications for galaxy clusters [7|, |8|, 
weakly magnetized non-rotating atmosphere local simulations have demonstrated that the 
MTI can ampli 
gradient |30l. 

& Stone [3l| studied the non-linear development of the MTI and HBI using numerical 



y magnetic field and lead to a substantial heat fiux down the temperature 
. McCourt et al. 24 1, Parrish & Stone 26 1, Parrish & Quataert j^^], Parrish 



simulations applied to clusters of galaxies. McCourt et al. |2J] found that the HBI can reduce 
the conductive heat fiux through the plasma due to the reorientation of the magnetic field 
lines. As the HBI should operate in the innermost ~100-200 kpc in the intracluster medium 
(ICM) of cool-core galaxy clusters, where the observed temperature increases outward, the 
cooling time of the ICM is shorter than its age; so the HBI removes thermal conduction as 
a source of energy for the cores, increasing the cooling fiow problem 28[. They also studied 
turbulence in the ICM and suggested that the interaction between turbulence and the HBI 
might be part of a feedback loop for the thermal evolution of the ICM [see also |29|, |38 |. 
These authors found that unlike the HBI, the MTI drives strong turbulence and operates as 
an effcient magnetic dynamo, much more akin to adiabatic convection, and that while the 
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MTI cannot saturate by reorienting the magnetic field, it can saturate by making the plasma 
isothermal. Sharma et al. {42 1 studied thermal instabilities including adiabatic cosmic rays. 
Balbus & Reynolds jsl showed that magnetic field configurations that nominally stabilize 
the HBI or the MTI can lead to g-mode overstabilities. 

Spherical accretion fiow is also subject to the MTI. In this sense Sharma et al. 43 1 
investigated the effects of MTI on spherical accretion flow using global simulations [see also 



111 ]. Consistent with previous local simulations, they found ampliflcation of the magnetic 



fleld and alignment of fleld lines with the radial direction (temperature gradient direction). 
Other scenarios such as the interiors and surface layers of white dwarfs and neutron stars 
have been also investigated [see e.g. Il2 |. 

Dilute magnetized plasmas not only exhibit anisotropic thermal conduction but also an 
anisotropic pressure tensor, which results from different kinetic temperatures for electron 
motion in the parallel and perpendicular direction to the magnetic fleld. For collisionless 
plasmas the pressure does not isotropize, and hence if we insist on using afluid description the 
standard magnetohydrodynamics theory (MHD) must be modifled 
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4Q|. 



The purpose of this work is to study the MTI and the HBI in a dilute magnetized plasma 
including the effect of pressure anisotropy. Our main goal is to determine the impact of 
these anisotropies on the instability criterion for the MTI and the HBI. To this end, we 
employ the theory of kinetic MHD (KMHD) as derived by the 16-moment method from the 
Vlasov equation p, |25|, |3^ . This formalism extends the well-known double-adiabatic theory 
of Chew, Goldenberg and Low (CGL) 13| in that it allows for nonvanishing (parallel and 
perpendicular) heat fluxes, which are neglected in the CGL theory. Moreover, the growth 
rates and instability criteria for hydromagnetic wave propagation obtained from the KMHD 
equations are in better agreement with those obtained from kinetic theory [see for instance, 



14 



16|. 



Although we keep the discussion quite general throughout the paper, we have in mind an 
astrophysical magnetized collisionless plasma, and as a typical example we consider the ICM 
(see section |V]). In particular, in the ICM one of the main open questions is the ampliflcation 
of the primordial magnetic fleld. This environment results unstable to the MTI on scales 
of ten kiloparsecs and larger outside cooling cores (regions where the temperature decreases 
outward). Parrish et al. 32| have shown that the MTI can produce convective motions 
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and a magnetic dynamo, leading to an efficient transport of heat. The analysis and results 
presented here may also be of relevance in similar astrophysical environments. 

The paper is organized as follows. In section |ll] we give a brief overview of the formalism 
we use to describe anisotropic plasmas. In section IIIII we describe the equilibrium state of 
the anisotropic plasma together with the linearized system of equations. We then present 
and discuss our results in sections IIVI and |Vl and finally we give a brief summary in section 

lYil 



II. BASIC EQUATIONS 



In this section we briefly review the formalism of KMHD and set the stage for the anal- 
ysis of wave instabilities given later on. More complete accounts of KMHD and similar 



approaches ( as we. 



in 16 
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1 as their relation to the kinetic theory of dilute plasmas) can be found 



40|, while an application of the 16-moment equations to the study of the 
firehose and mirror instabilities can be found in 

In a magnetized collisionless plasma such as those encountered in the interplanetary 
and intracluster mediums the pressure tensor is anisotropic with respect to the direction 
of the magnetic field. Thus, the transverse and longitudinal kinetic particle temperatures 
(associated with motion in the direction perpendicular and parallel to the magnetic field) 
will differ from each other, Ty 7^ T±. Due to the complexity of the kinetic equation for 
such a plasma, it is often convenient to employ a fluid description, which will naturally 
differ from the standard MHD description. In general, the single-fluid equations derived 
by Chew, Golderberg and Low (CGL) 13| are used. The energy equation of isotropic 
MHD is replaced by the double-adiabatic laws (or by the double-polytropic laws in the 
phenomenological approach of Hau & Sonnerup 19 | ]. Many studies of wave instability 
problems are based on the CGL equations [see e.g. [l8|, l20|. However, when compared to the 
results obtained from kinetic theory there is a discrepancy in the criterion for the slow-mode 
mirror instability, as well as basic differences between the nonlinear evolution of the mirror 
and firehose instabilities 

The inadequacies of the CGL approach stem from the unwarranted neglect of the third- 
order moment of Vlasov equation. The 16-moment method, which is a generalization of 
Grad's 13-moment method, provides a consistent way of deriving from the Vlasov equa- 
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tion the correct fluid equations for a heat-conducting anisotropic magnetized plasma. If 
the Larmor radius is much smaller than the other characteristic lengths of the plasma the 

25 1 [see also, [g], 



equations can be considerably simplified 



34 



. The equations describing 



the collisionless magnetized plasma then read 



1 + "^"' 







pg + 6V||Ap 
6(V-6) + V||6 



(1) 
(2) 



^+B{V ■v)-iB-V)v = 



where 



-V-Qii 



d 
dt 



(3) 
(4) 
(5) 

(6) 



Aj> = j>x b = B/B, V|| = 6 - V and Q\\,± are the parallel and perpendicular heat fluxes, 
which we assume for simplicity to be given by Braginskii's approximation [l^ 



Qll,± = -&XV||T|, 



(7) 



where x is the electron thermal diffusivity. Although x varies with temperature as 

X = 6x IQ-^T^/^ergs cm"^ K"^ (8) 



for simplicity and following Quataert 33| we will consider it as a constant in what follows. 
The specific entropies associated with parallel and perpendicular motion are given by [l| 

-In I '^'^ 

3 



and s 



^ln(^ 
3 \pB 



(9) 



where cy is the specific heat. Note that the total entropy 



s± = cv In 



1/3 2/3 

P\\ Pi 



p 



5/3 



(10) 



reduces to the ordinary MHD expression for the specific entropy in the isotropic case p\\ = p±. 
It will prove convenient for the problem at hand to employ the entropy production equations 
for sii ± instead of those for sii and s± separately. 



III. MTI/HBI IN THE PRESENCE OF PRESSURE ANISOTROPY 



In the following we ignore the ion contribution to the conductive heat flux, which is 



smaller than the electron contribution by a factor of ^Jmi/nie ~ 42. As stated in the 
previous section, we assume that the electrons have mean free paths much longer than their 
Larmor radius (as occurs, for example, in the ICM). Under these conditions the thermal 
conductivity of the plasma is strongly anisotropic. 

We consider a thermally stratifled plasma in the presence of gravity g = —gz, so that in 
equilibrium 

The magnetic fleld of the equilibrium state is assumed to be homogeneous and in the (x, z) 
plane 

B = B^x + BJ (12) 
so that there is a background heat flux given by 

dTw 

Q\\,± = -x{Kbz£ + blz)^ . (13) 



As in Quataert 



33|, we assume a steady-state initial equilibrium so that V ■ Q = 0, which 



implies a linear temperature proflle in z. Note that p± = p\\ + const, and that we consider 
a vanishing initial velocity. This equilibrium state is a solution to eqs. ([I])-(l5]). 

Under these circumstances, putting 6vz ~ exp^—iut + ik ■ r) and similarly for other 
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quantities the linearly perturbed equations can be written as 



uSp — pky = 



u[Kk - k'5v\ - {vj^ - Ap/p)k'b 



B 



-[kh - k%]h{5p±_ - 5p{) + i{kzk - k^z)g6p/p = 



u6B - Bk^ + {k ■ B)6v = 



■ Cv 
I — UJ 

3 



ds\\ 

+ 1 -r + a 

dz 



{5p\\ + 25p 
dsj_ 



dz 



5v,, 



.Cv 



1 



ds\\ 

dz dz 



{5p\\ - 25p, 
ds 



where we have used that 



6B 
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6p 
P 



(3 + 2a) 



~D,dlnT\\ dlnT^, 

■ ^XO—[—j — + a— J ) 

p dz dz 

Pf5T\\ ST A 

5B , ^ bp , 
— 1 + a + — 3-2a 
B p 



^D,d\nT» 

^XO—{—j — 
p dz 

P f6T\\ 

+X— \ -Tf- - 



rflnTi 



a- 



dz 



a- 



6J\ 
T, 



(14) 

(15) 
(16) 



(17) 



(18) 



SQi± = -x5&V||T|,,x - xH^b ■ VT^ 
va = B/s/^np is the Alfven speed, and we have defined 



(19) 
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IV. NEGLECTING PRESSURE PERTURBATIONS 



We will now consider the case in which pressure perturbations are neglected (Boussineq 
approximation). The growing modes of interest have growth times much longer than the 
sound crossing time of the perturbation, so it is sufficient to work in the Boussinesq approx- 



imation. In this subsection we closely follow Balbus |2| and Quataert 



33|. 



From eqs. (fT^ - (fT8|) in the Boussineq approximation we get the following dispersion 
relation 



where 



iAsu^ + A2UJ^ + {iAi + Ai)uj + Aq = (24) 



Ao = {u\ - u!)A2 (25) 

Ai = Ipcy^il + la){oo\ - F) - pN''^ (26) 

4 ~k± 

Ai = —gT\\Cv{l-a)b—{Kk,-b^k^) (27) 

A2 = p|j(l + a)w,,|| (28) 

A3 = ^pcyT|, (1 + '^a\ (29) 



and 



ua = k ■ va (30) 

= (31) 

^c,|| = Ix^b' (32) 
5 p\\ 

are the Alfven, Brunt-Vaisala and conduction frequencies, respectively, suitably modified 
by the anisotropy parameter, and F = (c^ y — c^^^)6^ with y ^ = p\\^±/ p. In eq. flMj) we 
have put 

i^' = ■1[(1 - 262)A;i + 26,6,A;,fcJ (33) 
for notational simplicity. To derive eq. ( IMl) we have used that Sp/p = —6T\\/T\\ = —STj_/T± 



see 



33|, and also that d\nT\\/dz = d\nT±/dz = d\nT/dz, where T = (T|| +2Tj_)/3 is the 
total temperature, which follows since a is assumed to be a constant. 

It is interesting to analyse the special case of a weak magnetic field, in which there exists 



an ordering in frequencies given by 



33] 



/ \ 1/2 

uJc':^ 0Jd= Ij^j > (34) 

where H is the local scale-height of the system and Ud its dynamical frequency. Besides, 
if a is sufficiently close to 1 then one can safely put \F\ < Uc- In this limit the dispersion 

8 



relation becomes 



~ gK — = gK— — (35) 

dz d z 



This expression is identical to the result of Quataert [33|, which shows that when the mag- 
netic field is sufficiently weak so that the timescale ordering given in eq. ( IMl) holds and 
a ~ 1 (as occurs e.g. in the ICM), the MTI and the HBI become independent of the 
pressure anisotropy. This is one of the main results of this work. 



Following Balbus [3[, Balbus & Reynolds [5| , we will now analyse the stability of solutions 
to eq. f l2^ using the Routh-Hurwitz criteria jl7|. It should be noted that the polynomial 
(Elj) has complex coefficients and therefore the generalized Routh-Hurwitz theorem must 



be used 4J]. Briefiy stated, the procedure to determine if a given complex polynomial f{z) 



with no roots in the imaginary axis is stable (Le. all of its roots have negative real parts) 



is as follows [for a detailed account see e.g. Il7|: (i) determine the real polynomials P{r 



and Q{r) such that f{ir) = P(r) + iQ{r) with real r; (ii) calculate the Sylvester matrix 
S{P,Q) associated to P(r) and Q{r); (ni) if at least one of the principal minors of S{P,Q) 
is negative or zero then f{z) is unstable. 

Applying this procedure to the polynomial in ( 12^ we find the stability criteria 



Aq > (36) 
^2 > and (37) 
^2(^1^2-^0) > Aj . (38) 



Comparing this result to the isotropic case studied in Balbus [3| , we see that the stability 
criteria get modified by pressure anisotropy in two ways. Not only do the expressions for 
the coefficients An contain a, but also an extra term Ai appears which vanishes if a = 1 or 
if b^kz = bzkx- The first effect does not change in a signiffcant way the quantitative analysis 
of Balbus 3|, so we shall not pursue it here. However, the appearance of Ai in the third 
stability criteria implies that when pressure anisotropy is taken into account A1A2 — Aq 
cannot be arbitrarily small, otherwise an instability develops. This represents a qualitative 
change with respect to the isotropic case, where A1A2 — Aq > for stable solutions. Note 
that this modiff cation only affects transverse perturbations, since Ai = if k± = 0. 
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V. APPLICATION TO THE ICM 



As an illustrative case, we will now present numerical results corresponding to the ICM. 
As typical values we take p ~ 10^^ cm'^, B ~ lO^^ G, T ~ 10^ K, and P ~ lO^^^ 
ergs cm~^. The value of g is estimated as the average of g{r) = GM/r^ from r = 
to r = 2rc, where G is the gravitational constant, M = is the mass of the ICM 

and Tc = 290 kpc. This gives (7 ~ 7 x 10^^ cm s~^. For the entropy gradient we take 
3d\ns/dz = —dlnT/dz, and taking the Brunt-Vaisala frequency to be ~ 10~^^ s~^ 



we get d\nT/dz ~ 10 cm . The anisotropy parameter in the ICM can be written as 



21 



37 



a = 1/(1 ±<y5), where (p = \p± —p\\\/p± ~ 1/vRe and Re is the Reynolds number 
To estimate the lower and upper bounds for the values of a we will consider Re > 50, so 
we get a G [0.87,1.16]. Note that in the ICM Re ^ 50 for convective motions, so that 
actually even lower values for |1 — q;| are expected in this environment. However, in order 
to show the physical behaviour of this instability as a function of the anisotropy parameter 
it is convenient to consider a broader, though unphysical, range for a (as we do in Figure 
[T]). In what follows we take a tangled magnetic field corresponding to = 0.3. 

In order to illustrate the dependence of the stability of solutions to ( 12^ on pressure 
anisotropy, we show in Figure [1] the growth time r = l/|Im(a;)| as a function of a for 
kz = kj_ = 1/kpc and = k± = 5/kpc. 

It is seen from Figure [1] that in the case of a = 1 the growth time, r, is of the order of 
~ 100 My 45|]. On the othe hand, the presence of anisotropy causes a strong decrease in r, 
which is considerably larger for smaller scales. The dependence of r with a is similar in both 
cases (a < 1 or a > 1). The inset shows the behaviour of r in the range 0.85 < a < 1.15, 
which as mentioned is a closer range to the values expected in the ICM. It is seen that 
for scales ~ 0.2 kpc the growth time in the anisotropic case can be almost one order of 
magnitude smaller than the one corresponding to a = 1. For larger scales (> 1 kpc) r 
decreases by a factor of ~ 2 or less. 

We now go over to analise the dependence of r with wave vector in the range [0.1, 10] x 
1/kpc (magnetic tension stabilizes the instabilities on shorter scales). Figure [2] shows r as a 
function of kj_ for kz = 1/kpc and for three values of a. We note that for other values of k^ 
in the studied range [0.1, 10] x 1/kpc the results are similar. It is seen that for any value of 
k± the growth time of the isotropic case is larger than the corresponding ones to a < 1 and 
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FIG. 1: (Color online) Growth time for the MTI/HBI as a function of a for = k± = 1/kpc and 
kz = k± = 5/kpc. The inset is a zoom for 0.85 < a < 1.15, which is a closer range to the values 
expected in the ICM. 

a > 1. The difference in r can reach one order of magnitude. As happened in the previous 
case, the behaviour of r is similar for both anisotropic cases. Interestingly, the decrease in 
r with the anisotropy is significantly smaller in the kpc range. 

We will now briefly discuss some possible implications of our findings to the ICM (in 
the case where the timescale ordering of eq. (|34l) does not apply). In this connection, the 
most important results of this work are that: (i) even with small pressure anisotropy the 
growth time of the MTI/HBI can become almost an order of magnitude smaller than the 
one corresponding to the isotropic case; (ii) the decrease of r with anisotropy is smaller 
for scales ~ kpc; and (iii) the effect of anisotropy on r is larger at smaller scales. As 
mentioned in the Introduction, the MTI/HBI can be a significant source of magnetic field 



amplification (ill . |24J . l31| . Our results imply that in the ICM the magnetic field amplification 
due to the MTI/HBI could be faster than what is expected in the isotropic case. Although 
we expect this effect to be small in the ICM because the magnetic field is weak and thus 
the timescale ordering holds, it could have some important consequences in the subsequent 
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FIG. 2: (Color online) Growth time for the MTI/HBI as a function of k± for kz = 1/kpc and 
a = 0.87; 1; 1.16. 

plasma dynamics following the linear regime. 

VI. CONCLUSIONS 

In this work we have studied the influence of pressure anisotropy, which is present in di- 
lute magnetized space plasmas, on the physics of the MTI and the HBI. We have performed 
a linear analysis based on fluid equations which go beyond the CGL double-adiabatic for- 
malism. 

The main conclusion that we can extract from our results is that, for the conditions 
prevailing in the ICM, the impact of pressure anisotropy on the MTI and HBI is small. 
More specifically, if the magnetic field is sufficiently weak so that the dynamical frequency is 
larger than any other frequency of the system, these instabilities do not depend on pressure 
anisotropy (provided the latter is small). 

On the other hand, if this timescale ordering does not apply the stability criteria for the 
MTI and the HBI will depend on pressure anisotropy in two ways, first through a dependence 
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of the terms which also appear in the isotropic case, and second through the appearance of 
new terms. We find that these extra terms affect the stabihty of the stratified plasma and 
the growth rate of the MTI/HBI. Specifically, we find that the growth time of the instability 
in the anisotropic case can be almost one order of magnitude smaller than the isotropic one. 

The analysis of the MTI/HBI in anisotropic plasmas presented in this work is largely 
idealizecL since we limited ourselves to the linear regime, we neglected viscosity and radiative 
cooling 5|, we considered static pressure anisotropy configurations (on which the MTI/HBI 
develops), and finally we used a simplified model of the ICM. In spite of this shortcomings, 
we believe that our analysis provides some insight, at least preliminary, into the effect 
of pressure anisotropy on the MTI/HBI, as well as on some possible implications for the 
dynamics of the ICM and similar astrophysical environments. 



Note added in v3 



Shortly after this work was submitted to the ArXiv, a preprint by M. Kunz 23| ap 



peared analyzing the impact of pressure anisotropy on the MTI/HBI by using the Bragin 



skii anisotropic viscosity equation for p_L — p\\ . In [23| , the author introduces time-scale 
and amplitude orderings measured by the Mach number M ^ 1, whereby the Boussineq 
approximation is implemented as Sp/p ~ 6T/T ~ 6p/{Mp) (i.e. relative changes in the 
pressure are much smaller than relative changes in the temperature or density), and not as 
6p± = 6p\\ ~ as we naively assumed here. 

This naive interpretation of the Boussineq approximation eliminates the Braginskii vis- 
cosity, and therefore in our work pressure anisotropy is not self-consistently calculated but 
rather acts as a (fixed) background on which the MTI/HBI develop. Due to this, our results 
differ from those of |23| in that we find, for the conditions prevailing in the ICM, only moder- 
ate changes in the MTI/HBI growth rates when including pressure anisotropy instead of the 



drastic changes found and described in [23[. This clearly highlights the necessity of including 
Braginskii viscosity self-consistently in order to study the effect of pressure anisotropy on 
the MTI/HBI as taking place in the ICM. 

In view of the above, we believe that the approximations involved in this work should 
be regarded merely as a first step towards a more thorough understanding of the MTI and 
HBI within the 16-moment formalism for collisionless magnetized plasmas. Concerning this. 
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the relation between the equation for the pressure anisotropy used in 23(] and the evolution 
equations for p± and of the 16-moment formalism is not clearly established. Hence it 
is still not clear what may be the consequences of the time-scale and amplitude orderings 
mentioned above on the 16-moment equations. We are currently investigating these issues. 
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